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Abstract 

The divided cell algorithm was introduced by Delone in 1947 to 
calculate the inhomogeneous minima of binary quadratic forms and 
developed further by E. S. Barnes and H. P. F. Swinnerton-Dyer in the 
1950s. We show how advances of the past fifty years in both symbolic 
computation and our understanding of homogeneous spectra can be 
combined to make divided cells more useful for organizing information 
about inhomogeneous approximation problems. A crucial part of our 
analysis relies on work of Jane Pitman, who related the divided cell 
algorithm to the regular continued fraction algorithm. In particular, 
the relation to continued fractions allows two divided cells for the same 
problem to be compared without stepping through the chain of divided 
cells connecting them. 

1 Preliminaries 

Notational conventions and a basic framework for working with approxima- 
tion problems are collected here for the convenience of the reader. 
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Diophantine Approximation problems deal with finding where the re- 
striction of a function to a special subset is small. In this paper, the function 
will be defined on the plane and the subset will be the integer lattice Z^. 
The function will take nonnegative values, so "small" will mean "close to 
zero" . It is customary to work with the reciprocal of the original function, 
and to freely treat oo = 1/0 as a number, since it will be a possible value of 
the supremum of a set of nonnegative numbers. 

By working with the name denoting a function instead of the traditional 
convention of denoting a sequence with subscripts, we are allowed the un- 
cluttered notation limsup/ to denote the infimum over all cofinite subsets 
5 of N of the supremum of the values of / restricted to S. 

The computation of the infimum of a function can often be organized 
by endowing the domain with a partial order for which the given function 
is order preserving. If the partial order has the property that descending 
sequences are finite, then the infimum of the values over the whole set is 
equal to the infimum over the set of minimal points for the partial order. 
This is valuable when the set of minimal points has special properties. In 
particular, it is usually possible to index the minimal points by the set of 
all integers Z so that a pair of adjacent minimal points has some special 
property. Such a function defined on set of all integers Z will be called a 
chain. 

The integer lattice in is 

identified with 7? by giving it a basis. Certain 
bases aid in the identification of the minimal points. These depend on the 
expression and have been called reduced. Dually, the expression giving the 
function in terms of a reduced basis has also been called reduced. 

Families of related problems lead to a space of reduced bases, and the 
study of all reduced bases for a single problem can be expressed in terms 
of a dynamical system on this space. This study will lead to strong results 
when the underlying space is compact. 

Our emphasis here will be visual with pictures of the plane including 
the lattice Z^. However, while a basis for the lattice is used in the algebraic 
description of the objects in the figure, other considerations may be used in 
the choice of viewing coordinates. 
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2 The MarkofF Spectrum 



Traditionally, as in [5], the homogeneous Markoff Spectrum is the set of 
values 

M(F)=sup| : x,yGZ,(x,y)/(0,0) | (1) 

where F{x,y) = Ax^ + Bxy + Cy^ is an indefinite binary quadratic form of 
discriminant D{F) = B^-iAC. The quantity M{F) is given as a normalized 
inverted minimum: normalized to allow natural comparison between values 
of different forms F; inverted to allow simpler expressions for interesting 
values in the spectrum. Those F with F{x, y) = for integers x and y (not 
both zero), as well as those taking arbitrarily small values, have M{F) = oo. 
The interesting cases are those for which M{F) is finite. 
Since F is indefinite, it can be factored over M. We write 

F{x,y) = {aQX + bQy){aix + biy), (2) 

and introduce new variables ^ = ooa; + h^y^ rj = aix + biy to get F = (^rj. 
Then, y^D{F) = \aobi - aibo\. 

In particular, the expression F{x, y) is encoded by the matrix 



ao bo 
ai bi 



(3) 



Left multiplication by this matrix takes the column with components {x, y) 
to one with components (^,7?)- Thus, it gives a change of variables between 
the arithmetic and geometric aspects of the study of the values of F on 
the integer lattice. The rows of the matrix are the coefficients in the fac- 
tors of F{x,y). As a change-of- variables matrix, its columns give the 
coordinates of the generators of the lattice. 

A change of basis in the lattice multiplies the matrix in ([3]) on the right by 
an integer matrix of determinant ±1; scaling the factors of -F(x, y) multiplies 
on the left by a real diagonal matrix. The value of M{F) is not changed by 
these actions. 

A visual approach to the Markoff Spectrum must show the integer lattice 
and the lines OjX + biy = = 0, 1). However, it is more convenient to use 
(.^, rf) as viewing coordinates since a fixed F will be studied using different 
bases for the integer lattice. In practice, this may be modified by a change 
of scale (^, rj) — > (a^, rj/a) in order to bring different lattice points into focus. 
Because of this choice of viewing coordinates, the lines where F{x,y) = will 
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Figure 1: Homogeneous example 



be called the axes of F. For example, a picture of F{x, y) = x"^ — 3y^ on the 
integer lattice uses viewing coordinates t]) with ^ = x + y\/^, t] = x — y\/3. 
Figure [T] shows this view of -F = (now just the coordinate axes), the lattice 
generated by {x,y) = (1,0) and {x,y) = (0,1), and the lattice cell whose 
(x, y) coordinates are (0, 0), (1, 0), (1, 1), (0, 1). 

In computing M{F), if a lattice point Pq is closer to both axes than the 
lattice point Pi is, then \ F\ is smaller at Pq than at Pi. This relation between 
Pq and Pi is a partial order of the type mentioned in the Preliminaries. 
Thus, only the minimal points for this partial order need be considered when 
finding M{F). Arranging the minimal points in order of their distance to 
a specified axis of F gives a chain of minimal points (except when an axis 
contains a nonzero lattice point). Since these results are well known, and 
have been given in detail in [3j, features of this chain are only sketched here. 
The corresponding results for inhomogeneous problems will be described 
later in more detail. Figure [1] shows that (x,y) = (1,0) and (x,y) = (1,1) 
are minimal points, but (x,y) = (0,1) is not since (1,0) is closer to both 
axes. 

A full description shows that two successive minimal points always gen- 
erate the lattice. These are the reduced bases of the lattice. In a precise 



definition of a reduced basis, it is convenient to fix the order of the axes, the 
order of the generators of the lattice, and to choose between a generating 
vector and its negative. With one set of choices, if x and y are the coordi- 
nates with respect to a reduced basis, one has ao > ai > and bi > —bo > 
in ([2|) and Since the matrix determines the reduced basis, it is ap- 
propriate to also speak of a reduced matrix when these conditions hold. The 
inhomogeneous case will also require matrices with 69 > > oq, but these 
are avoided in the tradition treatment of the homogeneous case. 

For three consecutive minimal points, a matrix whose columns are the 
basis consisting of the second and third points is the product of the corre- 
sponding matrix for the first and second points with the matrix 



1 

-a 



(4) 



with a = [oq /bo\ . To restore orientation and obtain the required signs of 
the matrix element, this must be multiplied on the left by a diagonal matrix 
whose first diagonal entry is positive and whose second entry is negative. In 
the example, when the matrix for the reduced basis (1, 1), (1, 0) multiplied 
by the matrix in dH) and rescaled, leads to the equation 



1 + ^3 1 ■ 




" 1 ■ 




i-Vs 1 




1 -2 





-1 + ^3 








2 

1-V3 



(5) 



The space of all reduced matrices with fixed determinant forms a com- 
pact set. Note that compactness requires that all inequalities be inclusive. 
Classical work often aimed for unique representations and required some 
inequalities to be strict, but sacrificing uniqueness to have a compact space 
of reduced matrices allows the Spectrum to be characterized in terms of 
attained extrema. 

The chain of the matrices given by ^ is one description of the steps in 
the continued fraction algorithm. It produces a symbolic dynamics that is 
useful for describing the relation between the reduced bases of a given form 
and the computation of M{F). In particular, a consistent choice of a vector 
from each reduced basis leads to a chain of minimal points {xn,yn), and 
M(F) = supM„(F) where 



M„(F) 



VD{F) 



(6) 



\F{Xn,yn)\ 

Each index n should be associated with both the minimal point (x^, y^) 
and the reduced basis with this point as first element. A sequence of in- 
dices can be found for which M„(F) M[F) and also the corresponding 
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reductions converge (see Lemma 6 of Chapter 1 of [5]). This shows that 
every value in the Markoff Spectrum is an attained supremum. This result 
is known as the Compactness Theorem for the Markoff Spectrum. 

A novel variation on this approach, allowing generalization to higher 
dimensions, can be found in [3]. 

The Divided Cell Algorithm transfers these properties of the continued 
fraction to inhomogeneous problems. 

3 The Lagrange Spectrum 

If the form F{x, y) in ([2]) is x{y — xa), then F(0, 1) = 0, giving M{F) = oo. 
However, if a is irrational, no other minimal points (x, y) have F{x^ y) = 0. 
If < a < 1, we set = (1,0) and (xo,yo) = (0)1)) giving a 

reduced basis; and then index the other minimal points by positive inte- 
gers. Properties of rational approximations to a are determined by L{q) = 
lim sup M„(i^) for n G N. The set of such values is called the Lagrange 
Spectrum. Theorem 1 of chapter 3 of [5] says that the Lagrange Spectrum is 
a subset of the Markoff Spectrum. This follows from the proof of the Com- 
pactness Theorem for the Markoff Spectrum. When L{a) is finite, the forms 
appearing in the construction are all equivalent to F, but the limiting form 
F* is nonzero on all lattice points other than the origin and L[a) = M{F*). 

In this paper, we will concentrate on the inhomogeneous Markoff Spec- 
trum, but applications to the inhomogeneous Lagrange Spectrum will follow 
by constructing a convergent sequence of reductions of the inhomogeneous 
expression x{y — xa — (3). 

4 The inhomogeneous Markoff Spectrum 

For inhomogeneous problems, the form F defined in ([2]) is replaced by 

Fi{x,y) = (aox + b^y + co)(aix + biy + ci), (7) 

while we continue to require {x,y) G 1? . Figures illustrating such problems 
will continue to be drawn in viewing coordinates for which Fi{x,y) = on 
the axes of the coordinate system. The origin of this coordinate system is 
no longer required to be a lattice point. For example. Figure [2] modifies 
the example of Figure [U by using factors = x + yV^ — 1 — 0.5\/3 and 
r]j = X — y\/3 — 1 -|- 0.5\/3 to study the expression Fj{x,y) = ^irjj. Note 
that the origin is now at (x,y) = (1,0.5). The parallelogram in the figure 
has vertices whose (x,y) coordinates are (0, 0), (1, 1), (2, 1), and (0,1) with 
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Figure 2: Inhomogeneous example 



edges that form a reduced basis. Earlier work has required that (cq, ci) not 
be in the lattice generated by (ao,ai) and (60,61) to explicitly exclude the 
homogeneous case. We propose to allow this case to exclude itself because 
it necessarily has a lattice point where Fj{x,y) = and interest will be 
centered on those Fj(x,y) that are bounded away from zero on the lattice. 
Notice that the parallelogram in Figure [2] has one vertex in each quadrant 
bounded by axes of Fj. This property defines a divided cell: the word "cell" 
refers to a fundamental parallelogram of the lattice, and it is "divided" by 
having its vertices separated by the axes. 

The definition of the inhomogeneous Markoff value is 



Mi{Fi) = sup 



«o6i - aibo 



Fi{x,y) 



x,y e 



(8) 



using the notation of d?]) . Note that the origin of the lattice is not excluded 
here, since it has no special role in Fj. 

To describe an inhomogeneous problem, the matrix of ([3|) must be re- 
placed by 

Oo 60 Co 
ai bi c\ 



B 



(9) 



Left multiplication by this matrix takes the column with components (x, y, 1) 
to one with components (C/i??/)- Again, the rows of the matrix are the 
coefficients in the factors of Fi{x,y). The interpretation of columns is a 
little different from the homogeneous case: the third column is the image of 
the origin, and the first two columns give generators of the lattice. If the 
matrix is augmented with a third row [0 1], one gets the affine change-of- 
variables matrix relating the column with components {x,y, 1) to one with 
components (.^/, 7]j, 1). In this matrix, a column with in the third position 
represents a direction; one with 1 in the third position represents a point. 
The vertices of the cell are found by adding the sum of a subset of the first 
two columns to the third column. All such columns have 1 in the third 
position, so they can be expected to represent points. 

Left multiplication by a two by two diagonal matrix changes the scale 
on the axes and right multiplication by an integer matrix with a third row 
[0 1] and determinant ±1 gives an afhne change of basis in the lattice. 

5 Divided cells as reduced objects 

We continue the convention of using the (x, y) for coordinates in the in- 
teger lattice, but describing geometric properties using (^/,r//) as viewing 
coordinates. 

Divided cells will be the reduced objects for the study of Fj{x, y) on the 
integer lattice. Several proofs of the existence of divided cells have been 
given, beginning with Delone [6] in 1947. Our proof uses work of Pitman 
[7], and will be given after discussing the role of divided cells in Diophantine 
Approximation. 

We choose the line aox + boy + cq = to be the vertical axis in Figure [2] 
with the positive halfspace on the right. Treating the vertices of the cell 
as the basic fundamental parallelogram of with (0, 0) in the lower left 
quadrant of the figure, gives cq < 0, 6o + cq < 0, ao + cq > 0, + 6o + co > 0. 
These inequalities imply oq > |6o|- A similar analysis of aix + biy + ci = 
as the horizontal line leads to 6i > |ai|. Conversely, these conditions on 
aQ,ai,bo, bi give a nonempty set of possible solutions for cq, ci. Since Mj{Ff) 
is invariant under scaling of the linear factors of Fj{x, y), one may introduce 
a convenient scaling, like oq = 6i = 1. Then ai and bo are each chosen from 
the interval [—1,1], and then each of cq and ci is chosen from an appropriate 
closed intervals. In this way, the space of divided cells can be represented 
by the fourth power of a closed interval. This scaling will not be used in 
this paper, but we will insist that ao > and bi > 0, forcing the base vertex 
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Figure 3: Cells and Boxes 



to be in the third quadrant. 

This construction shows that the specification of a divided cell can be 
done in two steps: first choose generators of the lattice giving the directions 
of the sides of the cell; then locate the origin. Barnes [1] introduced the term 
"I-reduced" for the lattice bases arising in this way. We keep the name, but 
take it to mean that oq > |ai| and bi > \bo\. 

If fli^o < 0, the cells are essentially the reduced cells of the homogeneous 
case. Such cells will be called Gaussian, or G-cells, indicating that they are 
reduced in the sense of Gauss. Note that, in contrast to the homogeneous 
case, no attempt is made to fix the sign of ai. 

The cells that are not Gaussian will be called N-cells. Since definitions 
should use inclusive inequalities, the correct characterization of an N-cell is 
dibo > 0. This allows a cell to be both a G-cell and an N-cell, but only when 
one of its sides is parallel to an axis. 

If a parallelogram is an I-reduced cell, then the possible locations of the 
origin in the cell form a rectangle inside the cell. This rectangle is called the 
inner box (which we will sometimes call simply a "box") of the I-reduced 
cell. Figure [3] shows two typical examples. In the figure, a G-cell is on the 
left and an N-cell is on the right. The width of the box is oq — so that it 
degenerates to a vertical line segment if ao = |ai|. Similarly, the inner box 
degenerates to a horizontal line segment if bi = \bo\. When both ao = |ai| 
and bi = |6o|, the box is only a single point. This illustrates that the first 
row, which gives the coefficients in the equation of the axis shown in the 
vertical position and describes the first coordinates of the cell, governs the 
divided cell step. 

As in the homogeneous case, the partial order defining minimal points 
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considers distances to both axes. However, this time it is necessary to treat 
each quadrant separately. Within a fixed quadrant, points that are closer 
to both axes will be smaller points in the partial order. We call this the 
basic partial order. It will need to be modified, but this is a good tentative 
definition. 

6 The divided cell algorithm 

Once one divided cell is available, it is possible to construct a chain of divided 
cells containing that cell. This construction is the Divided Cell Algorithm. It 
must be shown that Mj{Fj) can be computed using only the vertices of the 
cells obtained by this algorithm. This is essentially the content of Theorem 5 
of [2]. Another approach to using divided cells to compute Mj{Fj) is given by 
the theorem on page 530 of [6] . Our proof will distinguish six related chains 
arising from the divided cell algorithm: a chain of cells, a chain of boxes, and 
four chains of minimal points. The relations among these chains is not as 
direct as it is in the homogeneous case, so it is useful to keep them separate 
while showing how they are related. The chains of minimal points — one 
chain in each quadrant — play a key role in showing that all divided cells lie 
in a single chain and are used to characterize the quantity Mi{Fj) defined 
in dSI). Since the chains in different quadrants are independent, distances in 
different quadrants may be weighted differently. We don't explore that here, 
but some consequences can be found in Section 3 of [2]. Finally, the chain 
of boxes shows the simplest progression from one axis to the other. All of 
these chains terminate if there is a lattice direction parallel to an axis, but 
our statements will make no effort to distinguish that case. 

The construction of the Divided Cell Algorithm is used in two settings: 
given expression Fj, it produces its chain of divided cells; given only an I- 
reduced basis for a lattice, it describes all possible successor I-reduced bases. 

Theorem 6.1. If Fi admits one divided cell, then there is a chain of divided 
cells containing that cell. Given an I-reduced basis, there is one shape of an 
N-cell arising as a successor and one possible shape of a G-cell arising as a 
successor. The N-cell always occurs, but the G-cell may not: the number of 
positions of the N-cell is always one more than the number of positions of 
the G-cell. 

Proof. Suppose that we are given a divided cell defined by a matrix as in 
dl]). Since this is a divided cell, oq > |ai| and bi > I^qI, with additional 
bounds on each Cj in terms of and ftj. The details of the divided cell step 
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depend on the sign of 60. The algorithm terminates if 60 = 0, and there are 
only minor differences between the other cases, so only the case of 60 > 
will be illustrated. 

The definition of a divided cell then gives that cq < cq + bo < < 
cq + uq, so that the line segment from (co,ci) to (cq + 6o)Ci + 61) forms 
the left side of the cell and crosses the horizontal axis. This side can be 
extended until it crosses the vertical axis, giving an integer h > with 
Co + hbo < < Co + (/i + l)6o- The segment from T_ : (cq + hbo, c\ + hb\) to 
T_|_ : (co+/i6o+6o, ci+/i6i+6i) will form the top of the next cell. Similarly, the 
bottom of the next cell is found by extending the right side to get a segment 
from i?_ : (co+oq — A;&0i ci+ai — A;6i) to i?+ : (co+ao+6o^^^O) ci+ai6i — fc^i) 
for some /c > 0. Since segments r_T+ and B^B^ both cross the vertical 
axis, it follows that {h + k — l)bo <ao < (/i + /c + l)6o. This analysis shows 
that right multiplication by 



S 




1 





-1 

h + k 




1 



1 



(60 > 0) or 5 











1 

h+k 






1 - h 
1 



{bo < 0) 



(10) 

gives the matrix representing the next cell. In each case, h and k are positive 
integers with (h + k — l) \bo\ < ao < (h + k + l) \bo\. If 00/69 is not an integer, 
h + k must be one of the two integers nearest to |ao/6o|- 

The shape of the successor cell is determined by h + k; and the position 
by h. Those with different h and the same h + k are translates of one another 
The leftmost possible cell is an N-cell with h = 1 and k as large as possible. 
If |ao/6o| < 2, this is the only successor. Otherwise, the rules for determining 
the inner box show that decreasing A; by 1 and keeping h fixed gives a G-cell 
whose inner box abuts the box of this leftmost cell. Then, keeping this k 
and increasing /i by 1 gives a translate of the leftmost N-cell whose inner 
box abuts the box of this G-cell. The rightmost box will be an N-cell, and 
the the union of the inner boxes of these possible successors covers the inner 
box of the original cell. □ 

Note that the first row of ([9]), which gives the coefficients in the equation 
of the axis shown in the vertical position and describes the first coordinates 
of the cell, governs the divided cell step. 

The first part of Theorem 16. II is illustrated in Figure H] showing a divided 
cell with its box and, in two separate graphs, two successor cells with their 
boxes. In this picture, the original cell is a G-cell, and both types of successor 
are shown with the G-cell on the left (note that this figure contains G-cells 
with different signs of oi). Several lattice points are also included. 
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Figure 4: Successor Cells and Boxes 

Figure [H may also be used to analyze the chain of divided cell vertices 
in each quadrant. In the pictures, the common lattice direction of a cell 
and its successor gives a line joining the vertices of those cells in the second 
quadrant, and also in the fourth quadrant. Moreover, these two lines are 
adjacent lattice lines in that direction, so that there are no lattice points 
interior to the strip bounded by those lines. However, in the fourth quadrant 
of the second picture there is a point on one of these lines that is not a vertex 
of a divided cell although it meets our preliminary requirement for being a 
minimal point. We will now resolve this difficulty. We state the theorem 
for the first quadrant in order to have names for the edges that we use, and 
concentrate on points with small first coordinate but the proof is readily 
applied to both coordinates in all quadrants. 

Theorem 6.2. Given a divided cell C, let I be the projection of the open 
top edge of C on the horizontal axis. For each lattice line L parallel to this 
edge, let II he the points on L whose projection on the horizontal axis lies in 
I. Then each II contains at most one lattice point, and only those above the 
top edge of C contain such a lattice point in the first quadrant. Furthermore, 
the projections onto the vertical axis of the II are disjoint, so the ordering 
of these points by their second coordinate is the same as the order on the 
line L containing the point. 

Proof. For a lattice line L, the distance between consecutive lattice points 
on L is fixed, and the top edge of C gives one example of such a pair of 
consecutive lattice points. Again, since the lines are parallel, the difference 
of first coordinates is also fixed and equal to the width of / in this case. 
Hence, except when the endpoints of II are lattice points, there is a unique 
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lattice point in each 

Similarly, the relation between projections on the vertical axis of two 
consecutive II is also fixed, so it will be the same as the relation between the 
projections of the top and bottom edges of C. However, C is a divided cell, 
so all points on the top edge have positive second coordinate and all points 
on the bottom edge have negative second coordinate, so the projections of 
these edges are disjoint. For L below the top edge of C, all points of II have 
negative second coordinate, so II contains no point in the first quadrant. □ 

Any point in the first quadrant that is closer to the vertical axis than the 
vertex P of C in that quadrant must project into I, but Theorem 16.21 shows 
that all lattice points with that property have larger second coordinate than 
P. Hence P is a minimal point. 

When a side of C is extended to meet the positive vertical axis, one 
obtains a lattice line with one lattice point on each line parallel to the top 
edge of C. When the left side of the C is used in this construction, the 
first description of the divided cell step shows that the first lattice point 
in the first quadrant is a vertex of the successor divided cell. If it is the 
extension of the right side of C that meets the positive vertical axis, the 
first several lattice points will be in the first quadrant, but only the first 
and last of these are vertices of divided cells. This bypassing of minimal 
points in the divided cell algorithm is easily accommodated by augmenting 
the basic partial order. 

Theorem 6.3. If a line meets a quadrant in a bounded interval, the product 
of the distances to the axes is zero at the endpoints of the interval and has 
a unique interior maximum. The distance decreases as one moves from the 
location of the maximum towards either axis. 

Proof. A calculus exercise! When expressed in terms of one of the coordi- 
nates the distance is a quadratic polynomial with negative coefficient of the 
second degree term. □ 

When the line in Theorem 16.31 is a lattice line, this says that we may 
modify the basic partial order to also say that a lattice point is greater than 
another lattice point on the line that is on the same side of the point of 
maximum value of Fj and farther from that point. With this modification, 
the only minimal points on the line in this quadrant are the vertices of the 
original divided cell and its successor. 

Augmenting the basic partial order in this way on every lattice line and 
forming the transitive closure gives a new partial order called the extended 
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partial order with fewer minimal points. We will say the P is nearer that 
Q if P < Q m the extended partial order. 

Theorem 6.4. // divided cells exist, every minimal lattice point for the 
extended partial order is a vertex of a divided cell. 

Proof. This is now little more than using a known divided cell as the basis 
and using previous results of this section for an induction step. By symmetry, 
it suffices to show the result for a minimal lattice point P in the first quadrant 
that is closer to the vertical axis. By Theorem 16.21 there are finitely many 
minimal lattice points whose second coordinate lies between that of P and 
the original divided cell vertex in this quadrant. By the discussion following 
Theorem 16.31 the first of these is the vertex of a divided cell. There are 
fewer minimal lattice points in the first quadrant between this cell and the 
selected point, allowing induction to work. 

□ 

We illustrate the second part of Theorem 16.11 with Figure [5] showing the 
inner boxes of the possible successors. To draw both Figure H] and Figure [5l 
we used ao/ai = 2+\/5 4.236. The two parts of FigureOshow the number 
of each type of cell predicted by Theorem 16. 11 To avoid clutter, the cells are 
not shown in Figure [5l but the vertices are. Note that lattice points appear 
as vertices of the inner box of an N-cell, but the inner box of a G-cell is 
strictly interior to the cell and contains no lattice point. The cells that are 
collected in each of the pictures in Figure [5] are translates of one another in 
agreement with the expressions for their vertices appearing in the proof of 
Theorem 16.11 What Figure [5] shows is that the inner boxes of two successor 
cells intersect in at most a vertical segment and that the union of all of these 
boxes covers the inner box of the original cell. 

Combining this with the proof of Theorem 16.41 we find that the chain 
of boxes shows a systematic increase of height and decrease of width as we 
step through the chain. 

7 Pitman's Theorem 

Jane Pitman [7] related divided cells, which are the reduced cells of an 
inhomogeneous approximation problem, to the reduced bases of the corre- 
sponding homogeneous problem given by the continued fraction algorithm. 
Consequences of her work are an easy proof of the existence of divided cells 
(given here as the Corollary to Theorem lT.ip and tools for recognizing min- 
imal points. 
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Figure 5: All Successor Boxes 



Theorem 7.1. The cell of a Gaussian reduced form gives rise to two I- 
reduced N-cells. If the Gaussian cell has ai > 0, then the matrices corre- 
sponding to the other cells are obtained by multiplying its matrix by 
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The union of the inner boxes of these three cells is ( apart from duplication 
on the boundary) a fundamental domain for the lattice. 

Proof. Figure [6] gives a "proof without words" . It illustrates how the fun- 
damental domain of the given Gaussian cell may be cut into pieces that 
may be translated and reassembled to form the union of the three boxes 
described in the statement of the theorem. The parallelogram whose sides 
are not horizontal or vertical is the given cell. The box in the center of the 
figure is the inner box of this cell. The other boxes are the inner boxes of 
N-cells described in the statement of the theorem. The dashed line divides 
the part of the original cell outside the boxes into pieces congruent to the 
portions of the boxes outside the cell. □ 

Corollary 7.2. Every linear inhomogeneous problem in has divided cells. 

Proof. Employ the homogeneous theory to reduce the linear part (aox + 
b{iy){aix + hiy). Then locate the intersection of the axes in Figure [H The 
cell corresponding to that box is a desired divided cell. □ 
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Figure 6: Three boxes form a fundamental domain 

The N-cells described by Theorem 1 7 . 1 1 are called the neighbors of the G- 
cell in that theorem. One of these neighbors is characterized by |ao/6o| > 2; 
the other by |6i/ai| > 2. Conversely, each of these inequalities allows the 
construction of a neighboring G-cell of a given N-cell. 

If both inequalities hold, then the N-cell serves as an immediate link 
between consecutive reductions of the linear part of Fj{x,y). However, it 
is also possible to find N-cells for which neither of these will hold, so that 
they are not neighbors of a G-cell. Such cells will be considered in the next 
section. 

The boxes shown in Figure [6] give the matrices shown in (jlip . The third 
column affects only the location of a cell and not its shape, and is significant 
only for describing cells having some particular relation to the original G- 
cell. 

When ai < 0, the transition matrices of (jlip are replaced by 
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" 




1 





" 







1 





or 


-1 


1 





(12) 
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(If oi = 0, its sign should be chosen opposite to the sign of bo; if ai = 6o = 
the different constructions only involve cells with degenerate boxes.) 
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Figure 7: Superfluous cells 

8 Superfluous Cells 

This section investigates the role of the I-reduced cells that are neither G- 
cells nor neighbors of G-cells. We refer to such cells as superfluous cells for 
a reason that is given in Theorem 18. li 

Figure [7] shows the portion of a chain of divided cells starting with a 
cell C- for which —2 < ao/bo < —1 and bi/ai < —2 (the values used when 
drawing the figure were -(3 + \/5)/4 ^ -1.309016994 and -(3 + \/5) 
—5.236067977). The figure also includes the inner box of the first cell that 
is seen to also be the inner box of all cells shown. All but the last of these 
has a unique successor, and the figure shows this chain of unique successors. 
The last cell shown, C+, has ao/bo < —2, so there will be a choice of possible 
successors, none of which are shown. For all the cells C shown in Figure [71 
the cell C+ will be called the forward anchor of C and C_ will be called the 
backward anchor of C. 

Theorem 8.1. For a superfluous cell, the anchors are uniquely determined. 
For every vertex of a superfluous cell, a vertex of one of the anchors is nearer 
in the extended partial order. 

Proof. Since 1 < | ao/bo I ^ 2, the proof of Theorem 16.11 shows that the 
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divided cell algorithm involves a unique successor that is also an N-cell. As 
long as that cell is superfluous, the algorithm generates a unique forward 
chain. A closer examination of the function giving ao/6o for the successor 
in terms of the corresponding quantity in the original cell is an expansive 
mapping with ±1 as fixed points. From this, it follows that, apart from 
degenerate cases, the chain starting from any superfluous cell will reach 
a neighbor of a G-cell in a finite number of steps (this is not difficult to 
show, but the details are awkward to express, so they will be omitted). The 
process stops at the forward anchor of the superfluous cell. The process of 
stepping backwards through the chain of divided cells is governed by the 
ratio bi/ai in the same way, leading to the backward anchor of the original 
cell. Thus, Figure [7] describes the only way that superfluous cells can occur 
and relates these cells to the anchors. A study of the explicit matrix relating 
a superfluous cell to its successor shows that the anchors are attached to 
consecutive reduced bases of the lattice. 

Two of the vertices of a superfluous cell are also vertices of its inner box. 
These will be called the inner vertices of the cell. The inner vertices are 
shared with all cells shown in Figure [7] including the anchors, so they have 
now been found in a non-superfluous cell. 

The remaining vertices of the cells in Figure [7] (the outer vertices) lie on 
a lattice line parallel to and adjacent to the line joining the inner vertices, 
and one of the vertices of an anchor will be nearer in the extended partial 
order than a given outer vertex of a superfluous cell. □ 

9 A rigorous Framework 

The emphasis here has been visual. Figures were used to illustrate the con- 
structions and proofs. These figures were drawn using the Maple Symbolic 
Computation System. In order to tell the system what to draw, the cells 
and boxes were represented by matrices like S of Q. 

The visual approach was present in [6j, but was not used much by sub- 
sequent authors. Computers have facilitated the re-introduction of graphics 
into exposition, including the use of color where appropriate (the figures in 
this paper were presented in color at the conference). At the same time, 
increasing fluency in the language of Linear Algebra has encouraged the use 
of matrices to represent the objects met in the study. Our intent here was 
to use these developments to present old results in a way that will encourage 
new research. 

Some weaknesses of the Divided Cell Algorithm have appeared in our 
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exposition, but we have also shown that its appHcation to Inhomogeneous 
Diophantine Approximation can rely on methods like the ordinary contin- 
ued fraction that are associated to the Homogeneous Markoff Spectrum. 
Divided Cells become a tool for organizing the subject rather than a device 
for computing properties of individual problems. 
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